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Abstract 

The external-field method has been used extensively in the QCD sum-rule 
approach to explore various hadron static properties. In the traditional for- 
malism of this method, the transitions from the ground state hadron to ex- 
cited states are not exponentially suppressed relative to the ground state term 
and thus contaminate the ground state hadron property to be extracted. In 
this paper, we suggest a modified formalism, in which the transition terms 
are exponentially suppressed relative to the ground state term. As such, the 
pole plus continuum spectral model, traditionally invoked in QCD sum-rule 
approach, can be adopted. Thus, this modified formalism has potential to im- 
prove the predictability and reliability of external-field sum-rule calculations, 
which is illustrated in an explicit example. 
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The QCD sum-rule approach |IJ is a useful tool of extracting qualitative and quantitative 
information about hadron properties. One of the extensions of this approach made by Ioffe 
and Smilga |2| for external-field problems enables one to study hadron matrix elements of 
various currents and corresponding hadron static properties PHl8|. 

In the traditional external-field QCD sum-rule formalism, the contributions of transi- 
tions from the ground state hadron to excited states are not exponentially suppressed (after 
Borel transformation) relative to the ground state contribution which contains the ground 
state property to be determined. There are, in general, infinitely many transition terms as 
there are infinitely many excited states, and such terms should be retained in the calcula- 
tions. The usual approximation adopted in the literature is to introduce a new unknown 
phenomenological parameter (independent of Borel mass) accounting for the sum over the 
contributions from all the transitions between the ground state and the excited states. This 
new parameter is extracted from the sum rules, along with the ground state property. The 
isolation of the ground state property relies on that the polynomial (in Borel mass) behavior 
of the ground state signal is different from those of the transition terms. While this approx- 
imation has been used in earlier studies, it has been noticed recently that the parameter 
representing the transitions, in general, should be dependent on the Borel mass, which can 
have a sizable impact on the extracted ground state properties |TI|[r7|] . 

In this paper, we point out that a modified form for the linear response of hadron cor- 
relation function to the external field can lead to exponential suppression of the transition 
terms relative to the ground state term. As such, the pole plus continuum spectral model, 
usually adopted in QCD sum-rule approach, can be adopted. This formalism has poten- 
tial to improve the predictability and reliability of external-field sum-rule calculations. A 
comparison to the traditional formalism is illustrated in an explicit example. 

The external-field method in QCD sum rules starts from a correlation function of hadron 
interpolating fields in the presence of a constant external field |2|-|T8 



U(S,q) ee i J d 4 xe^(Q\T[ m {x)ri H m |0) s , (1) 

where S denotes the external field and t]h is a hadron interpolating field constructed from 
local QCD quark and gluon operators, carrying the quantum numbers of the hadron under 
investigation. The presence of the external field implies that 11(5', q) is evaluated with an 
additional term 



AC = -S ■ J , (2) 

added to the usual QCD Lagrangian, where J denotes the current of interest. Hereafter, all 
possible Lorentz indices of J and S will be suppressed. 

To proceed, one usually expands Eq. ([[]) to first order in the external field 

U(S,q)=U (q) + U 1 (q)-S + -- - , (3) 

where Ho{q) is the correlation function in the absence of the external field, which gives rise 
to the usual sum rules for hadron masses. One then focuses only on the linear response term 
Il 1 (g) which can be expressed as 
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U 1 (q) = d 4 xe^ x (0\T{rj H (x) 



d 4 yJ(y) 



%(0) |0) 



(4) 



The QCD representation of Ui(q), in terms of an operator product expansion (OPE), can 
be obtained by treating the external field perturbatively and then identifying the terms 
linear in S. The phenomenological representation for Hi(q) can be obtained by expanding 
Ili(g) in terms of physical hadron intermediate states or by invoking a dispersion relation. 
Here we follow the former to keep our discussion succinct. In a discrete-state approximation 
(i.e., neglecting the widths of hadron states, which is usually assumed in the QCD sum-rule 
approach [|IJ), the phenomenological representation can be expressed as 

Tli(q 2 ) ~ (0\ V \h ) (h \J\h ) (h \ri H \0) 



(M 2 



(0\ V \h ) (ho\ J\hi) (hi\fi B \0) +h.c 



{MS - ?)(M? - ?) 



+ terms involving only excited states 



(5) 



where \h ) and \hi) denote the ground and excited states with mass M Q and Mi, respectively. 
The first term contains the hadron matrix element, (ho\J\ho), of interest, the second term is 
the transition term from ground state to excited states usually discussed in the literature, 
and the rest involves only the excited states. After a usual Borel transformation one 
finds 



MM 2 )~§- 2 {h»\J\h Q )e 



-M 2 /M 2 



+ E 



\ \i(h \J\hi) 



Mi - Mi 



-M 2 /M 2 



exponentially suppressed terms 



(6) 



where Ao and Aj denote the coupling strengths of the hadron interpolating field to the ground 
and excited states, respectively, and M is the Borel mass. 

It can be seen clearly that the second term in Eq. is not exponentially damped 
as compared to the first term. This has been stressed constantly in the literature [Q-IS 



The usual strategy is to approximate the second (transition) term as Ae~ M ^ M2 with A a 
constant to be extracted from the sum rule along with the hadron matrix element of interest. 
However, one notices that A, in general, is a complicated function of the Borel mass M 2 . 
The neglect of such a Borel mass dependence will generate errors in the extracted hadron 
matrix element (see for example Refs. ]TU|,r7[). 



To improve this situation, we suggest to use the combination (Mq — q ) Ii\{q ), instead 
of Tli{q 2 ) alone. The corresponding sum rule for this combination has the form 



Xo\i{ho\J\hi) + h.c 



-Mf/M 2 



+ exponentially suppressed terms 

= \l(h \J\h )e^ M2 +5>e-^ 2 



(7) 
(8) 
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The transitions from ground state to excited states (second term) are now exponentially 
suppressed with respect to the ground state (first) term. This indicates that one may 
implement the traditional pole plus continuum model for the spectral ansatz of ITi. That is, 
one may model X^o e~ M ^ M2 in terms of the perturbative evaluation of III, starting from 
an effective threshold. As such, one no longer needs to introduce any phenomenological 
parameter to represent the transitions from the ground state to excited states as their 
contributions have been included in the continuum model. Therefore, this new formalism 
has potential to improve the predictability and reliability of the external-field method. 

As an example of illustration, we consider the proton matrix element of isovector-scalar 
current, H = (p\uu — dd\p)/2M p . This matrix element has been studied in Ref. |L6] within 
the traditional external-field QCD sum- rule method (see also Ref. Hl7| ). Here, our purpose 
is only to compare the new formalism suggested above with the traditional formalism. Thus, 
we shall adopt the same proton interpolating field and notations as used in Ref. [16] . 

In the traditional formalism, the sum rule used to extract H is given by [HJ 



p M 2 



-2aM 2 E L^ 9 + - X a 2 L 4 ' 9 + -mlaL 

3 A 3 3M 2 



-8/9 _ l fl 2j r 2/2r 



where M p is the proton mass, L = \n(M 2 / A 2 QCB ) / \n(fi 2 / A 2 QCD ) , and E = 1 - e - s °/ M2 



(9) 

with 
In 



so the continuum threshold. The values of a, tuq, Aqcd, and /i can be found in Ref. JTo 
the new formalism of Eq. (|8|), the sum rule becomes 

2HX 2 M p e- M H M " = -2a {M 2 p E - M 2 E,) M 2 L' AI% + - X a 2 M 2 ^ 

3 

+ \ml a M 2 p - | (l + ^ a 2 m 2 ^ , (10) 



where E\ = 1 — (l + s Q /M 2 )e~ s °/ M2 . In Eq. ([]), the transition term from the proton to excited 
states is not exponentially suppressed relative to the ground state term and is approximated 
as a constant phenomenological parameter A multiplied by e~ M vl M . This parameter is to 
be determined from the sum rule along with the matrix element H of interest. On the 
other hand, there is no such parameter in Eq. (|T0|) as the transition term is exponentially 
suppressed with respect to the ground state term [see Eq. (0)] and its contribution is included 
in the continuum model. 

To extract H from the above sum rules, we follow the numerical optimization procedure 

For defmiteness, we take for \ the value of x — 2.2 GeV -1 and the 
In the traditional method, the continuum threshold 



adopted in Refs. [16,17 



same Borel window as used in Ref. 16 



So is taken to be the same as that in the mass sum rules and H and A are extracted from 
the sum rule We use X 2 ~ 1.85 GeV 6 and s ^ 2.14 GeV 2 |L7||> which are slightly 



different from those used in Ref. |16| . In the new formalism, however, we extract both H 
and the continuum threshold So from the sum rule as the contribution of the transitions is 
effectively included in the continuum model. The predictions turn out to be 
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H ~ 0.60 [from Eq. (9)] , (11) 
H ~ 0.97 [from Eq. (0)] . (12) 

So, the two sum rules give very different predictions for the matrix element H. Which one 
is more reliable ? 

We have demonstrated in Eq. (^) that the parameter A in the traditional formalism 
should be a complicated function of the Borel mass. Ignoring this Borel mass dependence 
will alter the curvature of the phenomenological side of the sum rule and hence generate 
errors in sum-rule predictions. This, for the example under consideration, has been discussed 
in Ref. WJ\ (see also Ref. |T9|]). In particular, a complete form for the phenomenological side, 
which must be invoked in the traditional formalism, has been derived, and the corresponding 
sum rule takes the form |HJ 

2 H \ 2 p e' M p' M2 + A e- M p' M2 + B | e -°/^ 

= —2aM 2 E L~ 4 / 9 + *- X a?L^ + ^aL" 8 / 9 - . (13) 

It is the last term on the left-hand side which is ignored in the traditional formalism Eq. @. 
This term is Borel mass dependent and has impact on the prediction for H. There are three 
parameters H, A, and B to be extracted from this sum rule, and s is to be fixed at its 
value in the mass sum rules. Using the same numerical procedure, we find 

H ~ 0.98 [from Eq. (13)] , (14) 

which is very close to the result from the new formalism [Eq. (|i~2|)1. This indicates that 
the traditional treatment of the transition term has serious drawbacks arising from the 
approximation adopted, which can lead to large errors in the sum-rule predictions for the 
ground state matrix element. The new formalism, on the other hand, overcomes these 
drawbacks and greatly improves the reliability of the sum-rule predictions. 

One notices that the contribution of continuum in Eq. ( |i~0l) differs from that in Eq. @. 
This rises the question as to whether the change in the continuum contribution will affect the 
accuracy of the new formalism. We note that in the new formalism the contributions from 
the transitions are included in the continuum model. It is thus fair to make a comparison 
between the relative contribution of the continuum to the total phenomenological side in the 
new formalism and the relative contribution of the transition term plus continuum to the 
total phenomenological side in the traditional formalism. The result is shown in Fig. |l| We 
see that the contribution of the continuum in the new formalism is smaller than that of the 
continuum plus transitions in the traditional formalism in most of the Borel window. This 
can also be seen from that a term of the continuum contribution in the traditional formalism 
becomes two terms in the new formalism which tend to cancel each other [see Eq. (|]) and 
([10|)1. Thus, we conclude that the change in the continuum contribution does not affect the 
accuracy of the sum-rule prediction in the new formalism. 

Since the result from the new formalism is very similar to that from the traditional 
formalism with a complete form for the phenomenological representation [Eq. flI3"D], one 
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FIG. 1. The relative contribution (solid curve) of the transition term plus continuum to the 
total phenomenological side in the traditional formalism Eq. (||) compared to that (dotted curve) 
of the continuum to the total phenomenological side in the new formalism Eq. (|Iol). 



may wonder the advantages of the new formalism. An obvious one is that the usual pole 
plus continuum model can be used in the new formalism, which reduces the number of 
phenomenological parameters to be extracted from the sum rules [e.g., H and sq in the 
new formalism vs. H, A, and B in Eq. fll3|)1. A more important advantage is that the 
new formalism can be applied to external-field problems with any kind of external fields 
(vector, tensor, etc.) while a complete form for the phenomenological representation needed 
in obtaining Eq. (|H|) cannot be derived explicitly in the case of non-scalar external fields. 

The external-field method used in QCD sum-rule calculations concerns only the linear 
response term Hi{q) of the correlation function, which is actually a three-point function 
[Eq. (U)]. So one may evaluate Hi(q) by following the usual external-field method or by 
applying the OPE directly to the time-ordered operator product in Hi(q). The latter was 
initiated by Balitsky and Yung |2(J and followed by others [^T],^2j. The resulting sum rules 
in these two apparently different approaches must be identical. The apparent differences 
are only in the organization of the terms on the QCD sides of the sum rules. However, the 
calculations in the external-field method are usually simpler and more transparent than in 
the direct three-point function calculations. 

It is interesting to note that Balitsky and Yung [^UJ have proposed a procedure similar to 
the one discussed in the present paper. Such procedure has been adopted in later QCD sum- 
rule calculations based directly on three-point functions pl,22f, but not in the external- field 
method. However, the authors of Refs. pO| , pT| stated that the transitions from the ground 
state to the excited states are eliminated or canceled out in the sum rule for This 
statement is obviously not true. The transitions from the ground state to the excited states 
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do contribute to the sum rule, and their contributions are included in the continuum model. 



We also note that in Ref. [20-^] the continuum threshold in the sum rule for Hi(q) has 
been assumed to be the same as the one in the mass sum rules. This assumption, as shown 
above, is unjustified. Using the same threshold in these two cases will introduce artificial 
bias to the extracted hadron properties. 

As emphasized in Ref. [17|, it is the product of \ 2 H and the hadron matrix element of 



interest that usually appears in the sum rules for Ili(g). One then needs a good knowledge 
of A^- in order to extract the interested hadron matrix elements cleanly. This means that 
uncertainties associated with \ 2 H may give rise to additional uncertainties in the determina- 
tion of the hadron matrix elements. This is a general drawback of the sum rules for Hi(q). 
This drawback can only be sidestepped if the external field is treated non-perturbatively. 
The reader is referred to Ref. []23| for more discussions on this point. 

In summary, we have pointed out that a modified form for the linear response of hadron 
correlation function to the external field, concerned in the external-field method in QCD 
sum rules, can lead to exponential suppression of the transitions from ground state to excited 
states relative to the ground state term and hence reduce the transition contamination to the 
ground state property. This approach allows one to use the traditional pole plus continuum 
spectral model and thus has potential to improve the predictability and reliability of external- 
field sum-rule calculations, which has been illustrated in an explicit example. Our hope is 
that this new formalism of the external-field method will generate new interest in revealing 
hadron static properties via the external-field QCD sum-rule method. 
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